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1. Introduction 



Calculation of weak matrix elements of phenomenological interest is one of major applica- 
tion of lattice QCD. A calculation of four quark hadron matrix elements with the Wilson fermion 
encounters an obstacle since unwanted mixing is introduced through quantum correction with op- 
erators that have wrong chirality. 

One of the solution is to make use of the parity odd operator. By using discrete symmetries 
of the parity, the charge conjugation and flavor exchanging transformations it was shown ^ that 
the parity odd four quark operator has no extra mixing with wrong operators even without chiral 
symmetry. One of application of this virtue may be a calculation of the K Tin decay amplitude 
with the Wilson fermion. 

An improvement with the clover term is indivisible for the Wilson fermion. The RG improved 
gauge action of Iwasaki type has a good property at lattice spacing around a^' ~ 2 GeV imitating 
that in the continuum. It is plausible to use a combination of the Iwasaki gauge action and the 
improved Wilson fermion with clover term for our numerical simulation. Unfortunately renormal- 
ization factors of the = 1 four quark operators are not available for this combination of action 
except for AS = 2 part [^. A purpose of this report is to give the renormalization factor of four 
quark operators perturbatively which contribute to the K ^ nn decay. 



2. Four quark operators 

We adopt the Iwasaki gauge action and the improved Wilson fermion action with the clover 
term. The Feynman rules for this action is given in Ref. [^. We shall adopt the Feynman gauge 
and set the Wilson parameter r = 1 in the following. 

We shall evaluate the renormalization factor of the following ten operators 

eP«-l) = (^J)^ £ aj,"\qq)^,Q(^'') = (sxd), £ a^"\q X q)„ (n = 1,2,5), (2.1) 

q=u,d,s q=u,d,s 

Q(^"-') = isd\ X a^"\qq)^,Q^^"^ = (sxd), ^ a^^q X q)^, {n = 3,4), (2.2) 

q=u,d,s q=u,d,s 

a^'^ = (1,0,0), af = af = (1,1,1), a^^^ = a^^) = f 1,-1,-1) (2.3) 
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where 



{ld)j,n^ = sY^{\±Ys)d (2.4) 
and X means a following contraction of the color indices 

= X d)i^ {u X u)^ = {ladh)L {uhUa)L ■ (2-5) 

We are interested in the parity odd part only, which contribute to the K ^ nn decay amplitude 

^(2„-l) _ _ (2n-l)_ (2«-l) (2n) _ _Q{2n) _ M (n - \ 7 5) (2 6) 

^VA+AV ~ ^VA ^AV ' ^VA+AV ~ ^VA ^AV ' I"— K^-'^) 

„(2,i-l) _ „(2«-l) _ „(2n-l) „(2n) _ ^(2") _ ^(2n) / _ o (J n\ 

^VA-AV ~ ^VA ^AV ' ^VA-AV — ^VA ^AV ' I" — -'i^J, ' ) 
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Q^^^-'^ = {sd)y £ a(")fe)^, Qfy-'^ = {sd)^ £ «f'fe)y, (2.8) 

q=ii,d,s q=u,d,s 

QvA = X d)v L X ' Gav = X d)A £ a^'^ x q)y , (2.9) 

q=u,d,s q=u,d,s 

where current-current vertex means 

{ld)y {qq)j^ = {sYnd) {qYfiYsq) ■ (2-10) 

3. Renormalization factor in MS scheme 

(k) (k) 

We renormalized the lattice bare operators Q\J to get the renormaUzed operator We 
adopt the MS scheme with DRED or NDR. We notice there are two kinds of one loop corrections 
to the operators. One is given by gluon exchanging diagrams given in Ref. 10,01 for AS = 2 operator 
and the other is the penguin diagrams given in Ref. [|5|] for AS' = 1 operators. 

The renormalization of the operator is given by 

n^') _ yg n^J) I ypen^pen ysub^sub i n 

^MS~ O'^lat+'^i ^lat +A- P-^i 

where Q^^j^ is the four quark operators on the lattice, Q^^" is the QCD penguin operator and 01"^ is 
a lower dimensional operator to be subtracted. Zfj comes from gluon exchanging diagrams. Zf™ is 
contribution from the penguin diagram. 

3.1 Gluon exchanging diagrams 

For gluon exchanging diagram the one loop contributions are evaluated in terms of those to the 
quark bilinear operators by using the Fierz rearrangement and the charge conjugation Summing 
up contributions from three types of diagrams [Q] the one loop correction to the four quark operators 
is given in a form 

^ine-loop ~ -^7 Gtree; (3-2) 

where Qtril = Qva±av ^ level operator. The correction factors are already evaluated for the 
improved action in Ref. [Q] and is given as follows for our notation of the four quark operators 

T^lat T^lat T^lat T^lat T^lat T^lat 

''11 — ^22 —''33 —''44 —''99 —''10,10 

Tsf = T^f = l^(-^^H^cif + ^{Vv + VA)-:^(+Vs + Vp)], (3.4) 



„2 1 ^ 

^lat Tlat Tlat Tlat Tlat o /z:i„/"3 „\2 



= T^^ = T^^ = T^^ = T^^ = Tl^^^^ = -±-^-[6ln{Xar + Vv + Va-Vs-Vp) , (3.6) 



T^f = T^'f = ^^(-6ln{Xa)^-Vv-VA + Vs + Vp), (3.7) 
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where A is a gluon mass introduced for an infra red regularization and the number of color isN = 3. 
Vr is a finite part in one loop correction to the bilinear operator, which is evaluated in Ref. [^] for 
various gauge actions. 

The renormalization factor is given by taking a ratio of quantum corrections with that in the 
MS scheme multiplied with the quark wave function renormalization factor Z2 



Zfiina) = ^ ^ % ^ (3.8) 

Zfj{^a) = T^'-T^j (i^j). (3.9) 
The correction factor in the DRED MS scheme is given by 



T^MS 
''11 


T-MS 
— ''22 


T^MS 

_ — 


rpMS 7-MS 
I44 — 


T-MS 
— ''10,10 — 


7-MS 
■'12 


— ''21 


— ,^34 — 


T^MS T-MS 
''43 — ''9,10 


— ''10,9 — ■ 


T^MS 
''55 


T-MS 
— l-jj 


/N'-- 
~\ N 







y'^'^ (3.10) 



-3v^^ (3.11) 

(3.12) 

T^ = T^ = 3V^, (3.13) 

7^ = 7^8^ = 4^^^, (3.14) 

The same infra red regularization with the gluon mass should be adopted. The quark wave function 
renormalization factor Z2 is given in Ref. |^]. 
Substituting the above results we have 

Zfii^a) = Zf2 (Ata) =Zf3(Ata) =Z|4(Ata) =Zf9(Ata) =Zfo 10 (jua) 

Zl.i^a) = Z',,{^a) = \ + -l-^{-hn{^af+zl^ ' (3-17) 

Zli^ia) = Zl,{^a) = 1 + ^ (^^^^^^^ln(At«)'+4^ , (3-18) 

Zf^i^a) = Z|i(Aia) =Zf4(Aia) =Z|3(Ata) =Zf lo(Aia) =Zfo9(/ia) 

,2 

-31n(Ata)Vzf2), (3.19) 
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r,2 



Zf,(Ma) = Z%{^a) = (3ln{^ay + zl,) , (3.20) 

Zl,{na) = Zl,{^a) = ^zl, = 0. (3.21) 
The numerical value of the finite part is given in table [l| for A/^ = 3 as an expansion in csw 

^ij—^ii +(^SWZij +Cs^ffZlj ■ KJ-'^'^) 
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The finite part for the NDR scheme is given in table ^ We need to subtract the evanescent operators 
in the MS scheme, which comes from a difference of dimensionality from four for gamma matrices 
in operator vertex. 

Table 1: Finite part zf of the renormalization factor from gluon exchanging diagrams in the DRED scheme. 



^11 

(0) (1) (2) 


4 

(0) (1) (2) 


4 

(0) (1) (2) 


-23.596 3.119 2.268 


-25.183 5.420 2.923 


-18.041 -4.933 -0.020 



(0) (1) (2) 


7^ 

^56 

(0) (1) (2) 


-2.381 0.451 -2.020 


2.381 -0.451 2.020 



Table 2: Finite part z?' of the renormalization factor from gluon exchanging diagrams in the NDR scheme. 
csw dependent terms are the same as that in the DRED scheme (n — 1,2). 







^66 


^12 






-24.096 


-25.350 


-19.708 


-4.881 


-1.120 


-3 



3.2 Penguin diagrams 

Contribution from the penguin diagram is evaluated with the same procedure as in Ref . and 
the one loop correction to the four quark operators is given in a form 

/o(0 _ /-j,pen\lat pen fri 

tJone-loop ~ I i^tree> 

where is the penguin operator at tree level 

r,pen_/^(4) (6) \ _ }_ ( ,^(5) \ 

The correction factor is given by 

[TrT = 1^^ (ln«V + V^^l.) (3.25) 
with operator dependent factor 



C(!2i) 


= 0, c(e2) = i, 


C{Q3) 


= 2, 


(3.26) 




= C{Q6)= I 

q=u,d,s 


= Nf, 




(3.27) 


CiQs) 


= C(!27)=0, 




Nd 
2 ' 


(3.28) 


C{Qs) 


= c(eio)= I 

q=u,d,. 


OCq = Nu 


(3.29) 


C{Q9) 


= -1. 






(3.30) 
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p is a momentum of intermediate gluon propagator given in terms of external quark momentum, 
for which we set the on-shell condition. The finite part is expanded as 

,lat 



^pen = -1.7128 + CSW (-1.0878) . 
The correction factor in the MS scheme is given in a similar form 



n^'' _ /y^penyvi^ pen 

i^one-loop y i J i^tree' 



16712 3 

With the same infra red regulator p. The scheme dependent finite term is given by 



c(NDR)(2,) = ,(NDR)(g^^_^)^_l^ e(^DR)(22„)=0, {n>2) 



,(DRED) (g^) ^ ^(DRED) ^q^^_^^ ^ ^(DRED) ^q^^^ ^ 1 ^ > 2). 



(3.31) 

(3.32) 
(3.33) 

(3.34) 
(3.35) 



Combining these two contributions the renormalization factor for the penguin operator is given 



by 



pen 



npen MS ^ y,pen lat 



' pen 2 ' ' 



CjQi) 

16712 3 



lna2^2^^pen^^ 



(3.36) 
(3.37) 



Numerical value of the finite part is given in table I 

Table 3: Finite part of the renormalization ft 
0, 1) are given in the column marked as (k). 



Table 3: Finite part of the renormalization factor from the penguin diagram. Coefficients of the term c|jy {k = 



zP™(DRED)(«) 


zf"(NDR)(o) zP^,"_i(NDR)(0) zPf,"(NDR)(o) 


(,pen)(l) 


-0.2039 


1.0462 1.0462 0.0461 


1.0878 



3.3 Mixing with lower dimensional operator 

We shall evaluate the amputated quark bilinear vertex function given by 

,(sub) _ /„W 



%y = {Q'xM-p)dhp{p))^^^. (3.38) 

We consider a leading contribution to the vertex at tree level, which introduces mixing with lower 
dimensional operators. 
We immediately get 

itlvA = = «i"^5„.(75)„^ (/(-^)K) -/(-^)(m,)) , (3.39) 

= -Ci =N4"^5^,{rsU (/(-^)(m.) , (3.40) 
d^^l 4W{l,am) 



f (sub) / 



(3.41) 



(271)4 sin2/+W(/,am)2' 
W(/,am) =am + ^(l -cosZ^) . (3.42) 
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which may be evaluated with an expansion in the quark mass 

/(sub) 



1 d 



/(sub)(0) + 1^2l^^ 

d[am) a 2d[amy 



6 d{am) 



/(sub)(o) + ^(a) 



The numerical value is given by 



d{am 
d^ 



-/(™^)(0) = 
-/(™^)(0) 



1 



16712 
1 



;-21.466), 
(-14.92). 



(3.43) 

(3.44) 
(3.45) 



d{amY 167r2 

This contribution introduces a mixing with the lower dimensional bilinear operator ily^d) 
multiplied with a mass difference (m^ — nis). As it is clear from (3.42) this is due to the chiral 
symmetry breaking effect in the Wilson fermion. It may be better not to expand in quark mass 
since the coefficient ( 3.45 ) is rather large and d^ / d{amy' term has an infra red divergence 
at m = 0. The subtraction factor is given by 



4S = -2«f f/(™''('«^)-/' 



/(sub) I 



7 (sub) 

^(sub) 
■"2n-l 



-2Naf^ (;(sub)^^). 



/(sub) 



« = 3,4), 
(« = 3,4), 



^(sub) 
^2n 



:0, («= 1,2,5). 



(3.46) 
(3.47) 
(3.48) 



4. Conclusion 

In this report we have calculated the one-loop contributions for the renormalization factors of 
parity odd four-quark operators, which contribute to the K ^ nn decay amplitude, in the improved 
Wilson fermion with clover term and the Iwasaki gauge action. The operators are multiplicatively 
renormalizable without any mixing with wrong operators that have different chiral structures except 
for the lower dimensional operator. 
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